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bJ[). along one direction. As an illustration of the technique I construct cosmo- 



logical models which in their late time limit represent perturbations in the 



form of gravitational and scalar waves propagating on a non-static inhomo- 
geneous background. Several features of the obtained metrics are discussed, 
such as their early and late time limits, structure of singularities and physical 
interpretation. 
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I. INTRODUCTION 



The high degree of isotropy observed in the Universe on large scales today is usually com- 
bined with the Copernican principle to justify the assumption of homogeneity on the same 
scales. However, there is no known reason to assume that the Universe was isotropic nor 
homogeneous at very early epochs. The puzzling question of how the Universe might have 
evolved from an initially irregular state into the current isotropic and apparently homoge- 
neous state lacks a complete answer at present. To date, several regularization mechanisms 
have been put forward, such as Misner's chaotic cosmological program the standard 
inflationary scenario f||-|6| and, more recently, the alternative pre-Big Bang inflationary sce- 
nario @. However, none of these is completely satisfactory, and in general one cannot know 
for certain which range of initial conditions could have allowed the Universe evolve into its 
present form. Such scenarios can only provide one with partial indications of what initial 
conditions would have led a generic universe into the one observed at present. One way 
to maximize the amount of information obtainable from any cosmological program, such as 
those mentioned above, relies on studying the evolution of models with as many degrees of 
freedom as possible. This idea has motivated the special attention paid to inhomogeneous 
cosmological models in the last decades (see Krasinski || for a review). 

In general, attempts to obtain new inhomogeneous metrics involve some symmetry as- 
sumption, so that the full implications of the non-linearity of the theory are compromised to 
some extent. In addition, non- vacuum spacetimes are required to have a physically mean- 
ingful matter content, so that they portray realistic situations. Motivated by the possibility 
of the existence of non-trivial massless scalar fields in the early universe, I will concern 
myself here with cosmological solutions to the Einstein equations induced by such matter 
sources. In particular, I present a new algorithm to generate families of massless scalar 
field G\ cosmologies, i.e. time-dependent spacetimes with a single isometry. These new sets 
of metrics will be generated starting from generalized vacuum Einstein-Rosen spacetimes, 
which admit an Abelian group of isometries Gi acting transitively on spacelike surfaces. In 
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the last decades there has been intensive study of Gi vacuum and matter filled cosmological 
models and several major reviews on the subject have been written [P ,[T0 , 5[] . 



Given the large number of known G2 cosmologies and the various available techniques 
to generate further ones, algorithms transforming such spacetimes into G\ metrics represent 
powerful tools for generating new inhomogeneous solutions. The new algorithm, which I 
shall present, displays the nice property of reducing the symmetry while keeping the type 
of matter source unaltered. Nonetheless, the input and output solutions may not admit 
the same physical interpretation or even share some of the same relevant features. For this 
reason, if one wishes to grasp the physical meaning of every new solution, an independent 
analysis of it will have to be carried out. 

Here in I construct and analyze G\ models representing the propagation of gravitational 
and matter waves on a non-spatially flat background. The study of primordial inhomo- 
geneities in the form of waves is an active area of research. This is motivated by the fact 
that wave-like primordial perturbations originating from vacuum fluctuations during infla- 
tion may be responsible for structure formation. Unlike other types of inhomogeneities 
formed in the early universe, they would have remained nearly unaltered up to the present, 
and therefore allowing the possibility of their detection. 

The exact inhomogeneous G\ spacetimes seeded by gravitational and matter waves stud- 
ied here represent a generalization of more symmetric configurations considered by Charach 



and Malin [11]; in those models the background hosting the waves was homogeneous. In the 



context of colliding plane waves, G2 diagonal spacetimes with waves of scalar and gravita- 



tional nature have also been considered. Space-times such as those studied by Wu [12 



or 



any solution generated by the methods of Barrow JT3 and Wainwright et al. [14j] could be 



taken as starting points to construct new G\ metrics modeling interactions of waves on a 
curved background. 

Furthermore, interest in this generation procedure is not restricted to the relativistic 
framework; solutions to Einstein's equations with a massless scalar field (in what follows 
m.s.f.) may be used to generate solutions to alternative theories of gravity, such as Brans- 



Dicke theory or string theory in its low energy limit. In the latter case, one could even take 
those spacetimes to generate new solutions with other massless modes in the characteristic 
spectrum of the theory. 

The plot of the paper is as follows: First I introduce the G\ massless scalar field solution 
generating algorithm itself. Then, I construct new inhomogeneous metrics starting from a 
infinite dimensional family of solutions which in the WKB limit admit an interpretation in 
terms of waves. It will be shown that at early times these solutions behave like a Belinskii- 



Khalatnikov generalized model |L5[ with homogeneity broken along one spatial direction. 
The structure of spacelike singularities of the new solutions in the early time limit will be 
analyzed as well, and the special features due to the presence of the matter source will 
be indicated. Next I will consider the solutions' high frequency limit and show that they 
can be thought of in the same physical terms as their G 2 counterparts. In particular the 
new cosmological models represent a spacetime with a spatially inhomogeneous background 
filled with a non-static inhomogeneous scalar field and a null fluid of "gravitons" and "scalar 
particles." As time grows the null fluid's contribution to the energy-momentum tensor grows 
faster than that associated with the homogeneous part of the scalar field generating the 
background geometry. Thus, the model evolves into an inhomogeneous generalization of the 
cosmological model of Doroshkevich, Zeldovich and Novikov (DZN) ||16|| . Finally, the main 
conclusions are outlined. 



II. SOLUTION GENERATION ALGORITHM 

Basically, the new generating technique is a prescription to break homogeneity along one 
direction in G% m.s.f. cosmologies, ending up with new spacetimes possessing a single isom- 
etry but having the same type of matter content. Remarkably, the pioneering investigations 
on matter filled universes of Einstein-Rosen type, carried out respectively by Tabensky and 
Taub [|17]] and Liang [|11J considered models with non interacting scalar fields, even though 
at the time there was no clear physical motivation. 
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According to a conjecture due to Belinskii, Lifshitz and Khalatnikov (BLK) fl^-24 



G 2 metrics seem to be specially relevant for the description of the early universe, as such 
solutions give the leading approximation to a general solution near the singularity at t — 0. 
Their claim has very recently found the support of numerical results ||25| . In particular BLK 
considered approximate Einstein-Rosen solutions and performed an analysis of their local 
behavior in the early and late time regime. An interesting result reached in the course of 
those investigations, which is specially relevant for the present paper, was the prediction of 
a high frequency gravitational wave regime in the late epochs of the Universe. 

Before going any further it is convenient to explain how G2 spacetimes induced by a 
m.s.f. can be generated starting from vacuum solutions to Einstein's equations with the same 
symmetry. For the sake of simplicity, the discussion will be restricted here to a particular case 



of a well-known general procedure p6|,p!]]J^,p!3|J14,28,29|. The generic diagonal line-element 



with G2 symmetry will be taken starting point: 

ds 2 v = e Mt ' z) (-dt 2 + dz 2 ) + G v (t, z) (f^ z) dx 2 + e~ p ^ z) dy 2 ) , (1) 

where the subscript v stands for vacuum. A new solution g^ u of the Einstein equations with 
a massless scalar field <p as a source and line-element 

ds 2 = e f{t ' z \-dt 2 + dz 2 ) + Git, z) (e p ^ z) dx 2 + e'^dy 2 ) , (2) 

can be obtained by the following transformations: 

G = G V , (3a) 

p = Bp v + C \ogG v , (3b) 

f = f v + Ep v + F\ogG v , (3c) 

<p = Ap v + D \ogG v ; (3d) 

provided that the constants A, B, C, D, E and F are subject to the constraints: 

BC + 2AD = E, (4a) 



C 2 + 2 D 2 = 2 F , (4b) 
B 2 + 2A 2 = 1. (4c) 

The conditions (|j) arise by demanding the following are satisfied: 

Rfiu = WVV > ( 5a ) 
^V M V^ = . (5b) 

In principle a large number of new m.s.f. G<i cosmologies can be obtained by simply 
applying the procedure sketched above to any of the representatives of the populated family 
of vacuum Einstein-Rosen spacetimes. However, generating m.s.f. solutions with a lower 
degree of symmetry is a more cumbersome task. At this point I would like to draw attention 



to a method given by Feinstein et al. [|30| which allows one to generate families of solutions 
with a two-dimensional degree of inhomogeneity and a self-interaction term for the massless 
field of the form 

V = V (\)e- x *. (6) 

With these procedure new metrics are obtained by introducing an x-dependent conformal 
factor on the input G 2 metric, and where in general the potential term only vanishes for 
|A| =6. This difficulty in canceling the self- interaction term for the scalar field can be traced 
back to the high degree of symmetry of the x-dependence in the models considered in ref. 
0- 

In order to find a prescription for introducing an additional degree of symmetry in m.s.f. 
Gi metrics without switching on a potential in the process, I have considered the possibility 
of having a more general x-dependence in the metric. In particular, I have sought metrics 
of the form: 

ds 2 = Q(x) e nt ' z) (-dt 2 + dz 2 ) + G(t, z) (e p{t ' z) dx 2 + E(x) e- p{t ' z) dy 2 ) , (7) 
and made the following ansatz for the scalar field: 



<p(t,z,x) = (p(t,z)+A(x). (8) 

Note that the massless scalar field case included in the solutions of Feinstein et al. is also a 
particular case of the models here. Since the requirement is that no potential should arise 
in the transformations, the equations that must hold are: 

Riw = $,\jp,v , (9a) 

^V M V^ = 0. (9b) 

Explicitly, equation (j9a|) is equivalent to the set of equations: 

#oo = #oo + e J p = (p jt , (10a) 

Ru = Ru - e J p = <p tt , (10b) 

#oi = #oi = <P,t<P,z , (10c) 

#22 = #22 + ^ — + — = V,x , (10d) 

#12 = — ' ' H L r= L - = ip z ip x , ( lOe) 

R ° 2 = Yen + ^T = ' (10f) 

#33 = #33 + e p — e p = ; (lOg) 

whereas equation ([TFJ) in explicit form reads 

Inspection of the equations (|10| ) indicates that the three ^-dependent metric functions must 
have the form: 

n(x) = x k , (12a) 

E{x) = x n , (12b) 

A (a;) = mlog \x\ ; (12c) 

subject to the following constraints on the parameters k, n and m: 



k(2k + n-2) = 0. 
n(2k + n-2) = 
k + nC = 2Dm, 
nB = 2 Am, 
2 m 2 = 4 k - 3 k 2 . 



(13a) 
(13b) 
(13c) 
(13d) 
(13e) 



In this case, the Klein-Gordon equation ([5b|) reduces to 



m (2k + n-2) 







(14) 



which is automatically satisfied provided that ( |13a| , |13bj ) and ( |13ej ) hold. Note that the 
parameter k must be non- negative and not larger than 4/3. The ^-dependent term of the 
scalar field will be maximum for k = 2/3. Consistency of the solutions is reflected by the 
fact that for the following choice: 



m = n = k = 



(15) 



the set of equations (jT^JTJ) is satisfied for any value of A, B, C and D, and the input G 2 
m.s.f is recovered. On the one hand, for m ^ and k 1, if k and C are taken as free 
parameters, one can parameterize the solution's constants in the form: 



n = 2 - 2k. 



m = sign m 



A = v^sign (A) 



V4k-3k 2 



V2 
1 - k 



B = sign (A m) 



2-k 
1-k 



2-k 



V4:k-3k 2 



1-k 



D = sign (m 



k + (2-2k)C 
V8 k - 6 k^~ 



E = sign (A m) 



2-2k 



2-k 



C (2 - k) 2 + (2 - 2 k) k 



y/4k-3k 2 



F _C 2 (k + (2-2k)Cf 
2 8k -6k 2 



(16a) 
(16b) 

(16c) 

(16d) 

(16e) 

(16f) 

(16g) 



Four different subcases can be distinguished, depending on the choice of the sign of A and 
m. On the other hand, in the particular case m ^ 0, k — 1, the constants take the values: 

V2\m\ = y/2\D\ = \B\ (17a) 

A = (17b) 

C 2 + 1=2F (17c) 

E = sign(B) (17d) 

In general the metric admits only one Killing and one homothetic vector, namely: 

d 

Zke = ^ (18a) 

thv = Vj-, (18b) 

where by definition: 

C-Smfa = 0, (l 9a ) 

£tto9w = 2 9^ (l 9b ) 

It is known that a matter source and the geometry it induces need not share the same 



symmetry properties unless it is a perfect fluid |31| , |32| . This property is usually referred to 
as the "inheritance problem". Inspection of the metrics obtained here show that for k = 4/3 
the x-dependent term of the scalar field is absent, even though the metric depends on that 
coordinate in a non-trivial way. This is equivalent to: 

C imc T^ = (20a) 
A™ C V ± 0; (20b) 

where 

OX 



and the vector £ mc is a so-called matter collineation [10,33,32 



Bearing in mind the similarity between the new algorithm and the one of []3U| , one might 
wonder whether geometries like (|7|) can be also seeded by a scalar field with an exponential 



potential. Such a situation may be shown to be only possible if n = k, which is nothing 
but the case already found by Feinstein et al. In order to prove this, let us consider the 
case where the generic geometry ([?[), under the constraint ([12]), is induced by an exponential 
potential V{(p) = V Q (X) e~ x ^. In this case the field equations are: 

Kv = + 9„vV e~ X( P , (21a) 

V'V^-^M (21b) 
dip 

It is only necessary to look at the equations for i? o an d -R33 to realize the following constraint 
must hold: 

k(2-2k-n) =n(2-2k-n) = 4V ^ 0. (22) 

Compatibility of the latter set of equations in the case of a non vanishing potential requires 
n = k, or in other words, that the ^-dependence of the metric is given by a global conformal 



factor as in the case studied by []30 1 . 



It is important to note here that the generation technique does not restrict the character 
of the gradient of metric function G(t, z) of the input metric. That function determines the 
local behavior of the spacetime, and its gradient can be globally timelike, spacelike, null or 
vary from point to point. Although in what follows I am focusing on a case with a timelike 
character of the gradient of G(t, z), a window is left open for the study of other physically 
appealing cases. 

Moreover, even though our generating prescription has been used to break homogeneity 
of an input m.s.f. solution with a single degree of inhomogeneity, it is also possible to 
construct an equivalent algorithm that would transform certain static spacetimes into non- 
static ones. One should start with a m.s.f. model with two commuting Killing vectors, one 
of them being timelike, and then generalize it by introducing time dependent factors in the 
metric and scalar field as I have done here. 
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III. COSMOLOGIES WITH GRAVITATIONAL AND SCALAR WAVES 



After having outlined our method to generate uniparametric families of G\ cosmologies, 
I shall now construct the counterparts of a family of G2 cosmologies found by Charach and 



Malin [11|. That set of metrics, which can be thought of as inhomogeneous sinusoidal pertur- 
bations of the well-known Bianchi I spacetimes, represent the propagation of gravitational 
and scalar waves on an anisotropically expanding flat background. Simpler models of that 
sort, not including scalar degrees of freedom, were studied earlier on by Berger P5| . 

In this paper I am only considering Charach and Malin's solutions in the asymptotic 
limits t ~ and jt 3> 1 for any value of j. It is convenient, however, to give here the full 
expressions of the metric functions of the vacuum solution from which they were derived, 
namely: 

G v = t , (23a) 

00 

p v = plogt + cosfj'O - Zn)]ajJ (J t) , (23b) 
i=i 

fv = — - — + flJ2 a 3 C0S \. a i( z ~ z j)] J o(jt) + — f{[ajJo(jt)] 2 + [ajJ^jt)} 2 } 

z j=i z 1=1 j=i 1 j 

x {sin[l(z - zi)]sm\j(z - Zj)] [laciaij Ji(lt) J (J t) - jaiacjJ Q (lt) Ji(jt)] 

+ cos[l(z - zi)} cos[j(z - Zj)] [jaiajJi(lt)J (jt) - laiajJ (lt)Ji(Jt)]} . (23c) 

Charach and Malin's solutions were obtained by breaking the spatial homogeneity of 
Belinskii-Khalatnikov homogeneous solution ||15|| , which also has a m.s.f. as a seed. The 
gravitational and scalar degrees of freedom of those spacetimes satisfy linear wave equations, 
with the form of cylindrically symmetric waves propagating on Minkowski spacetime, for 
which the spatial and temporal coordinates have been interchanged. Since the solutions are 
of the standing wave form they are fully compatible with the 5* 1 ® S 1 ® S 1 topology (three- 
torus). However, the G\ counterparts of those spacetimes cannot have the same topology, 
x cannot be a cyclic coordinate in this case due to the presence of the term proportional to 
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log \x\ in the scalar field (p. 



A. Early time behavior and singularities 

From the analysis of the new G± solution's early time behavior it can be determined 
whether spacelike singularities arise at the time origin t = 0. A look at expressions ( |23|) 
shows that in the case under discussion, the periodic inhomogeneities can be neglected in the 
very first stages of that spacetime's evolution. Our metrics will be then G2 inhomogeneous 
generalizations of the cosmological models of Doroshkevich-Zeldovich-Novikov (DZN) ||16|| , 
which have three commuting Killing vectors. In this limit the metric and scalar field read 

~ diag(-x fc ei(t),x fc ei(t),e 2 (t),x n e 3 (t)) (24a) 
(p = ipo \ogt + mlog \x\ (24b) 



where 



ei = t fo , (25a) 
e 2 = t 1+P0 , (25b) 
e 3 = t 1 -* , (25c) 



and 



Po = BP + C , (26a) 
f = ^^ + E(3 + F, (26b) 
y?o = A(3 + D . (26c) 



In addition, the following relationship holds: 

fo = ^±+A (27) 



Following Charach and Malin, the metric can be rewritten using a synchronous set of coor- 
dinates in which the new time coordinate r is defined 
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d T = y/ei(t)dt; (28) 
this way the metric can be recognized as a simple inhomogeneous generalization of a 



Belinskii-Khalatnikov |15| solution. 

In broad terms, the presence of a spacelike singularity at early times will be reflected in 
the behavior of the curvature invariants R, R^R^ and R^ ua iR^ ua& . Due to the inhomo- 
geneous character of the metric it is possible in principle to have a conspiracy between the 
parameters, so that on certain hypersurfaces some those invariants are identically null and 
therefore do not reveal the presence of a singularity in the spacetime. Since the Kretchmann 
scalar K = R^ U(T sR' lua6 is however non-null and positive everywhere I will use this in the 
search for singularities. In the case here it explicitly reads: 



J£ _ 3/q 1"/- . OnO . „ / „ . .lOl . > 



t6/0+ 4 x6fc lv 1 + Po) 2 + 2 (/o + iy\ + 16t2/0+4x2fc [(i - P oY(i +f+ P y 

+(l +p f(l + f-pf} . (29) 



Since (p7) holds it can be seen that K is singular at t = for any value of the three free 



parameters of the solution, indicating thus the generic presence of a spacelike singularity at 
the time origin. 

It is also possible to study the singularity structure of the solutions in a more refined way, 
in particular by studying the expansion along each spatial axis. In general, the behavior will 
strongly depend on the values of the parameters k, C and (3. It can be shown analytically 
that for large enough (3 there will necessarily be contraction along the z axis, moreover if 
k > 1 that will be the case regardless the value of j3 and C. Another fact that can be easily 
checked is the impossibility of having simultaneous contraction along axes x and y. Three 
main types of singular behavior can be distinguished: 

a) Point-like singularities (Quasi-Friedmann behavior) 
All three spatial directions shrink as the initial time t = is approached; or explicitly 
lirrit->o €j = Vi. Depending on how many directions have the same expansion rate the 
behavior will be completely anisotropic, axially symmetric or isotropic. 
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b) Finite lines 

This type of singular behavior occurs when in the vicinity of t = one of the spatial directions 
neither expands nor contracts with time; in other words, it is said that the direction i is a 
finite line if Iim t _,o£j = 1- The subcases can be classified according to which the direction 
behaves in that way. In general, there will be a single finite line, though it is possible to 
have particular cases in which a second finite line exists. However, it is impossible to have 
such behavior along all three directions. 

c) Infinite lines (Quasi Kasner regime) 

An infinite line along the i direction exists when lim t ->o Q = oo. Again, three cases can be 
seen to occur, depending on which is the axis displaying that feature. For some particular 
values of the parameters the maximum allowed number of two infinite lines can be reached. 

Since the singular behavior of the metrics depends on three parameters it becomes rather 
complicated to represent graphically the structure of singularities in the general case. For 
that reason I will restrict myself to cases in which the value of one of the parameters is fixed, 
namely C. In the pictures here two different sets of lines can be distinguished. On the one 
hand, the black lines in each plot correspond to the curves along which a spatial direction 
becomes a finite line. In the region delimited by the black continuous line and the axes, a 
infinite line type of singularity arises along direction z. Here one can see how for C = 0.5 it 
is possible to have simultaneously the same behavior along directions y and z. For k > 1 in 
the region above the black dashed line there is contraction along direction y, this behavior 
gets reversed for k < 1. Similarly, for (5 values less than those along the dashed-dotted 
line contraction takes place, and the contrary happens for A; > 1. The points where two 
black lines intersect correspond to having two finite lines. On the other hand, the grey lines 
represent the curves along which two spatial directions display the same expansion rate. 
The fact that at a given point the three grey lines intersect reflects the possibility of having 
isotropic expansion, and for the two C values chosen, that point lies in the k > 1 region. 
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FIG. 1. Structure of singularities of the inhomogeneous generalization of Belin- 
skii-Khalatnikov's model for C = 0.5 (l.h.s. figure) and C = —0.5 (r.h.s. figure) with sign (Am). 
On the one hand, the black lines indicate the value of (3 as a function of k for which a given direc- 
tion behaves as a finite line. The black continuous, dashed and dashed-dotted lines correspond to 
a finite line along z, y or x directions respectively. On the other hand, the grey lines indicate the 
value of [3 as a function of k for which two spatial directions have the same expansion rate. The 
grey continuous, dashed and dashed-dotted lines correspond respectively to e 2 = £3, e\ = 63 and 
ei = e 2 . 

B. WKB limit 

With regard to the physical interpretation of the new solutions constructed here, it is 
their high-frequency limit which turns out to be most interesting. This limit, also called the 
WKB regime, corresponds to the regime in which the time elapsed since the beginning of 
the Universe is much larger than the period of any perturbation mode. By taking jt ^> 1, 
for every value of j, in the normal mode expansions of the scalar field and metric functions, 
Charach and Malin were able to show that the relativistic solutions taken here as input, 
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represent scalar and gravitational waves propagating on an spatially flat background. In 
this limit, such universes are causally connected because the particle horizon is larger than 
the wavelength of any of the modes of the independent degrees of freedom, namely the 
transverse part of the gravitational field p and the scalar field (p. 

In this vein, it will be proved here in two different ways that the G\ counterparts to 
Charach and Malin's cosmological models can also be thought of in terms of waves propa- 
gating on a non-static background. Thus, the physical interpretation is not spoiled in the 
process of homogeneity breaking. A direct consequence of the additional degree of inhomo- 
geneity present in the G\ solutions, as compared to their more symmetric counterparts, is 
that the background on which the waves live is not spatially flat. Cosmological backgrounds 
perturbed by waves have also been considered in a series of paper by Centrella and Matzner 
|3~iH36| who studied collisions of plane gravitational waves in these settings. 

Let us consider now the late time expressions for the metric functions and the scalar field 
of the Gi spacetimes obtained by applying the new technique to the cosmological models 
given by expressions (|2"B|): 



p ~ Pologt + Bp(t, z) , (30a) 
/~/ log*+(^— J£jaj, (30b) 
(p ~ y2 logt + Ap + mlog |x| , (30c) 



00 



2 7T 

P = ^\j ~~ t a i cos ^ 1 ~ 4) C08 \j( z ~ • ( 30d ) 

Since in this regime p < 1, the metric can be split into a background fj^ plus a 
perturbation metric h^, that is: 

9(tv ~ fjnv + , (31a) 

77^ = diag (-x k t fo e flt ,x k t fo e flt ,t 1+Po ,x n t 1 - po ) , (31b) 

= diag (0,0, t 1+p0 Bp, -xH^Bp) . (31c) 
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Here, in addition to (0), I have made the following definition 

/i = E^(^ + 2^ 2 ). (32) 

3=1 271 

A peculiarity regarding the perturbations on the scalar and gravitational degrees of 
freedom is that for < k < 1 they will be on phase whereas for 1 < k < 4/3 they will be 
phase-shifted by tt. Nonetheless, it will be seen later that whatever the value of k the scalar 
and gravitational perturbations contribute constructively to the energy momentum tensor. 

A proof of the assertion made above regarding the non spatial flatness of the spacetimes 
described by the background metric fj^, is provided by the expression of the spatial curvature 
of the t = constant three dimensional hypersurfaces: 

2 

(3,7?) p_ 171 

Let us now proceed to analyze the tensor and thereby show that it represents the 
gauge-invariant perturbations of a background spacetime fj^; or in other words, that it 
satisfies the wave equation V^V^h^ = 0, or equivalently, following f3"7|| , that: 

h!y = , (34) 
VU 7A = , (35) 

with 7, A = 2, 3, and where the D'Alambertian and the covariant derivative must be calcu- 
lated using the corresponding background metric f]^. 

Since it is straightforward to see that the trace-free condition fl3~4]) is satisfied by con- 



struction, the problem reduces to satisfying (35), which is equivalent to requiring: 



1 Though the factor B 2 + 2 A 2 equates to unity in the simple case I am dealing with, it has 
been deliberately introduced in the definition of /i ; so that the trail of the separate contributions 
to the energy momentum tensor of the graviton and scalar field pair can be followed. Had I 
considered the general case of the procedure to generate a G2 massless scalar field solution, then 
/i/E^iJ« J V(2vr). 
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h 2 2,x - 2 T\ 2 h 2 2 = , (36a) 

^-2^33 = 0. (36b) 

On the one hand, expression ( |36a| ) is identically null because neither 7/ 22 nor /i 22 are 

^-dependent . On the other hand, it can be seen that ( |36 b| ) is also satisfied by just hav- 
ing in mind that: 

h 33 = B f) 33 , (37a) 

fla = flog Jiks) . (37b) 



Once it has been proved that the wave equation V 7 V 7 /i M ^ = holds one can properly refer 
to the tensor h^ u as describing metric perturbations in the form of gravitational waves. 
In order to give additional arguments in favor of the interpretation of the solution in 



terms of waves propagating on a non flat background, I shall follow Charach and Malin fll 
and analyze the energy- momentum tensor of the background metric fj^. It will be shown 
that the stress-energy tensor is naturally manifested in terms of two components. One of 
these corresponds to a null fluid, supporting thus the interpretation suggested above; while 
the other term corresponds to an inhomogeneous massless scalar field with no ^-dependence. 
In particular, 

{ V )fu = (1) fv + (2) fu ^ ( 3g ) 

where (1) f^ ^ and (2) f; = 0. Explicitly 

Wf ° = -r^^- - t -(2+fo) 1 + 2 fo - Po 2 (39a) 



= -t- {1+po) h + r(2+/0) 1 7e"7° ' (39b) 

«2 ? = r a +A )2a^zlh± j (39c) 



(2 + /o) 1 + 2/0-P0 2 



(i)f 3 = f -(i+po) m , f -(2+.h) 1 + 2 fo - PO 

3 2 a; 2 4e/i'x fc ' 



(39e) 
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( 2 )T ° = -r( 1+ ^^V-, (39f) 
(2)j-i = r d+/o)_^ . (3 9g ) 

I will now proceed to give an interpretation for the ^'T* term. The Klein-Gordon 
equation for a scalar field ip, calculated using the metric fj^, takes the form 



X 



- 0, (40) 

l-k efltffo-po V ' 



a solution of which is 

■0 = v?o logt + m log |x|. (41) 

The energy momentum tensor for the field ip propagating on the spacetime fj^ yields 
exactly. That being so, one can conclude that the exact solution obtained after applying 
the generating technique to ( P3"D asymptotically evolves into a solution with a single degree 
of inhomogeneity; that is, the sinusoidal inhomogeneities along the z-axis vanish with time. 
In other words, ^'T" corresponds to the energy-momentum tensor of the exact solution one 
would obtain by switching off the periodic inhomogeneities along z. 

On the other hand, the traceless term ^T" can be shown to account for waves. It can 
also be separated into two parts, namely: 

(2) T F ^ w 'f F+ Wf Fl (42a) 
iGW %» = t (42b) 

j=— oo \J I 
oo J^2 a 2 

^ — X] 2\j\ K>1 i Ku i' (42c) 

The null vector is defined by 

K w - = -i=(U|,j,0,0). (43) 

V7TI 



It is clear that ^T^ u corresponds to a null fluid describing a collisionless flow of "gravi- 
tons" and "scalar particles" . It is easy to see that in the case under discussion gravitational 
waves will be absent only if the metric has no x-dependence. 
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I has been shown that the new solutions represent, in their WKB limit, waves propagating 
on an inhomogeneous non-static spacetime. In our case the only difference with respect to 
the case studied by Charach and Malin is that the background is not spatially flat. 

As far as the further evolution of the model is concerned, the presence of the additional 
degree of inhomogeneity in the model plays a crucial role. The null fluid's contribution to 
the energy-momentum tensor dominates the one due to the homogeneous part of the scalar 
field, and in the t >> 1 limit the energy-momentum tensor has a term accounting for the 
waves, plus another corresponding to the ^-dependent term in the scalar field. The scalar 
curvature of the background in this limit will be given by: 

2 

Mr — (44) 

so that this universe can either become singular at t = oo or flat, depending on the sign of 
1 + pq. The possibility of having a curvature spacelike singularity at late times is entirely 
due to the inhomogeneous character of the background. 

Another proof of the interpretation of the solutions in terms of a null fluid propagating 
on the background fj^ can be given. Let us calculate the energy- momentum tensor that 
corresponds to the scalar field (p in the late time limit and just retain terms up to the order 
t -1 . Under this restriction the only non-null terms of the energy momentum tensor are Too 
and Tn, which are given by: 

/ A 2 \ °° °° /— ( TT TT 

Too = - f n = y—J J2Ys a i a iV l j |sin[K* - z t )} sin[j(z - Zj)\ cos[Zt - -] cos[jt - -] 
-cos[l(z - zi)] cos\j(z - Zj)] sin[/t - J wa\jt- J j + C(t" 2 ). (45) 
Averaging T 00 over a region 0<t<2 7r, 0<2<27ron the (t, z)-plane one obtains: 

i 

3=1 

In the same approximation it can be seen that 



(Too) = -(T n > = — - 2 / / (T 00 )dzdt = — -^jA 2 (46) 

4 TT JO JO lilt ~\ 



(T) = <T|> = (T|) = 0(r 2 ). (47) 
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So, essentially it has been found that 

(T^> =^ f, u . (48) 

Summarizing, both methods have yielded the result that for t ^> 1, and up to the order £ _1 , 
the energy-momentum of the m.s.f. can be reduced to the null fluid form. 

The WKB regime of the solutions under discussion admits a reformulation in terms of the 
density of particles contributing to the modes of two fields. I shall strictly follow Charach's 
approach here P8] , which consists of performing a quasi-classical treatment based on the 



geometrical optics energy-momentum tensor. A family of Lorentz local frames is introduced 
so that the density of particles in each normal mode can be defined through: 



= \i\Hw\ ffl i ( no summation over v) (49) 

where fj uu represents the components of the inhomogeneous generalization of the DZN metric. 
A set of observers corresponding to this tetrad are characterized by the 4-velocity: 

u v = Afa = ( Jff°, 0, 0, 0). (50) 



Let us consider now equations (42^, 42c), which give the WKB stress-energy tensors of the 



cosmological model with gravitational and scalar wave perturbations, namely: 



B 2 a 2 
A 2 a 2 



( GW )j' lxv j — 3 K^jKvj (51a) 



( sw )T /iv j _ ^ 3 K^jK U j. (51b) 
The density of scalar and gravitational particles in the n-th mode is given by 

f(0)GW 2 2 

where k^ v is a null vector with dimensions of length. Besides, since the description here 
is based on units c = G = 1, the Planck constant has dimensions of (length) 2 , where 
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h ~ 10 66 cm 2 . The density does not depend on the direction along which the particles 
propagate, that is: 



In the G\ models considered here, the volume of the spatial sections t = constant is not 
finite; for that reason the total number of particles in each mode of the two degrees of 
freedom has no upper bound. 

In light of this reformulation it was suggested that one should regard the evolution of G 2 
models filled with waves as describing a process of transforming the initial inhomogeneities 
along z into quanta of various fields. Clearly, this interpretation's validity is extendible to 
models with just one isometry, such as those constructed in this paper. 



Before I finish, I will summarize the main results. I have presented the first method to 
generate uniparametric families of general relativistic cosmologies having a two-dimensional 
inhomogeneity and a m.s.f. as a source. In the context of either General Relativity or 
alternative theories of gravity, one can generate a large number of new inhomogeneous 
cosmologies using this algorithm, where moreover the only input needed is any of the many 
known vacuum relativistic cosmologies with two commuting Killing vectors. 

It has also been shown that this technique allows one to construct families of cosmologies 
which represent waves propagating on a spatially curved cosmological background. The 
spacelike singularity structure of these solutions has been studied, and several peculiarities 
due to the matter content have been elucidated. 




(53b) 




IV. CONCLUSIONS 
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